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Abstract 

A physical mechanism has been proposed to describe the formation of a stable state 
of a nonwetting liquid after filling of a porous medium at an increased pressure with the 
subsequent reduction of the overpressure to zero. It has been shown that the observed 
transition of the system of nonwetting-liquid nanoclusters to the stable state in a narrow 
range of filling factors and a narrow temperature interval is due to a decrease in the surface 
free energy without change in the chemical potential owing to the collective interactions 
between liquid nanoclusters in the neighboring pores. This effectively corresponds to their 
mutual attraction. The model makes it possible to describe the observed dependences of 
the volume of the liquid in pores in the stable state on the filling factor and temperature. 

To interpret the data on the characterization of the porous medium by means of mercury 
porosimetry [IrH], it is necessary to understand the mechanisms of phenomena and processes 
occurring when a porous medium is filled with a nonwetting liquid. One of these phenomena 
is entrapment, i.e., the formation of a stable state of the nonwetting liquid filling the medium 
at the increased pressure with a subsequent reduction of the overpressure to zero. Information 
on the distribution of the confined dispersed liquid over pores of various sizes is also necessary 
for, e.g., the interpretation of investigations of the dependence of the melting temperature on 
the size of clusters in the case of entrapment [HJE]. The stability of the dispersed nonwetting 
liquid in the porous medium was observed in a number of works studying the intrusion of 
mercury into porous glasses [HE], as well as the intrusion of water and aqueous solutions of 
salts [71412] and organic materials |4^ [T3TfT6] into various hydrophobic silica gels, structured 
porous media MCM-41, and zeolites. It was found that the fraction of the volume of a liquid 
confined in a medium can vary from to 1 in media with different porosities and depends on 
the average radius, width of the pore size distribution, concentration of substances in water, 
and temperature. As far as we know, the stability of a dispersed nonwetting liquid in a porous 
medium has not yet been systematically studied. For example [UHZMlH] that this phenomenon 
can be due to the structural features of the medium: the pore size distribution, the barrier 
for the formation of the liquid-vapor interface at extrusion from a filled cylindrical pore, and 
hysteresis of phenomenological wetting angles at intrusion-extrusion. Analysis indicates that 
all mentioned data cannot be described under these assumptions. 

In [19], it was shown that the fraction of a liquid in a stable state depends on the filling 
factor and the transition of a part of the liquid to the stable state occurs in a narrow range of 
the filling factor and in a narrow temperature range. 

In this work, taking into account both the percolation properties of the medium with pores 
of various sizes and the "multiparticle interaction" between liquid clusters in neighboring inter- 
connected pores, we propose a mechanism of the appearance of a stable state of a nonwetting 
liquid in a porous medium. The energy of the "interaction" between two liquid clusters in 
neighboring pores connected through a mouth (throat) is defined as a change in the surface 
energy of the liquid in the mouth after the appearance (disappearance) of a liquid cluster in the 
neighboring pore. The stable state of the nonwetting liquid appears as a result of a dispersion 
transition when filling factors and temperatures are such that the energy of the "multiparticle 
interaction" (attraction) between neighboring liquid clusters is higher than the energy of the 
liquid-solid interface. The transition to the stable state is described in the framework of the 
analytical percolation theory. 



The fraction 62 of the volume of water remaining in the L23C8 porous medium after the 
first intrusion cycle as obtained in |19| is shown in Fig. [I] as a function of (a) the fraction 9\ of 
the volume of pores filled in the first intrusion-extrusion cycle for three temperatures 279, 286 
and 293 K and (b) the temperature at d\ = 1.0. 
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Figure 1: (a) Dependences 62 (#1) for three temperatures T = 279, 286, and 293 K and (b) 
dependence #2 CO a t @i = 1-0 (experimental points in comparison with the calculations using 
the model proposed in this work shown by the lines) 

According to Fig. [U the fraction of the volume of the remaining liquid at 9\ > 0.16 is 
02 ~ 0.06 at T = 293 K. A decrease in the temperature to 286 K does not change the remaining 
volume at Q\ < 0.85, but leads to an increase in the fraction of the confined liquid at 81 > 0.85 
(62 ~ 0.6 at 61 = 1). With a further reduction of the temperature to 279 K, the fraction of 
the confined volume remains unchanged at 9\ < 0.85 and a further increase in the fraction of 
the remaining liquid is observed at d\ > 0.85. Thus, the dependence of the fraction of the 
volume of the remaining liquid on both the fraction 6\ of the volume of the pores filled by 
the liquid and the temperature of the system has a threshold behavior. It follows from Fig. [T] 
that the threshold values of the temperature and fraction of the filled volume are T c w 287 K 
and 6 c q « 0.85, respectively. In this work, we propose a model that takes into account both 
the properties of the interface and the percolation properties of the space of the pores in the 
medium and makes it possible to describe the observed phenomenon. 

We assume that a porous body is formed by a solid framework in which pores form a spatial 
structure as overlapping spheres distributed in the radius. This model of the porous medium 
is a generalization of the model of randomly located spheres, which is widely used to describe 
porous media |20j . In the model of randomly located spheres, an elementary pore is represented 
as a spherical cavity with notches (throats). In the model under consideration, a "quantum" of 
a change in the volume of the liquid in the medium at quasistatic intrusion (extrusion) is the 
volume of one pore. It is assumed that the volume of the throats is negligibly small compared to 
the volume of the pores. When the liquid is intruded into a pore and is extruded from it, liquid 
menisci are formed it these throats. We suppose that the width of the pore size distribution 
5R satisfies the inequality 5R/R < 3, which guarantees that the spread of the percolation of 
the threshold in interconnected pores is independent of the radius of pores [21J. 

Pores can be filled only when they are connected with each other through throats and with 
the surface of the porous medium. This occurs if porosity is such that the system of pores in 
it is above the percolation threshold, i.e., ip > ip c , where ip c = 0.16 -j- 0.3 is the percolation 
threshold [20| 122] . which characterizes the porous medium. Pores are connected with each 
other because an infinite (geometric) cluster, which consists of interconnected pores, appears 
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at ip = (p c . At (p > (p c , the intrusion of the liquid into the disordered porous medium is the 
intrusion of the liquid into the percolation cluster consisting of interconnected pores of various 
sizes. Only pores with the radius larger than a certain value can be filled at a given pressure. 

The extrusion of the liquid from a pore in such a medium with a decrease in the overpressure 
can be considered as a result of the extrusion of the liquid from a pore belonging to the 
percolation cluster of interconnected pores of various sizes filled with the liquid. The liquid can 
be extruded from the pore under two conditions. First, its extrusion should be energetically 
profitable; this means that negative work should be spent on the extrusion of the liquid from the 
pore owing to a change in the energy of the porous medium-liquid interface and the formation 
of menisci in the throats of the neighboring pores. Second, since the liquid can be extruded only 
through the connected system of filled pores, geometrical paths should exist for its extrusion 
from a given pore. 

Entrapment of the liquid in such a medium can be due to two reasons: (i) geometrical, when 
clusters containing a finite number of filled pores without paths for extrusion are formed in the 
porous medium and (ii) energetic, when the liquid in pores remains in the stable state with a 
decrease in the overpressure to zero because of the positive work necessary for the extrusion of 
the liquid. Figure [2] shows the schematic representation of a change in the state of the porous 
medium and the number of menisci in the pore before and after extrusion at various filling 
factors. Figures |2K I2t, and [2fe show the states before the extrusion of the liquid from the pore 
and Figs. [2b, |2H, and[2f show the states after extrusion of the liquid from the pore. Black color 
indicates filled pores, white color marks empty pores, and shading indicates the framework of 
the porous medium. Arrows in Figs. [2hr-[2]i indicate menisci and their numbers z and k. In 
Figs. |2fe and [2JT, 1 is the "infinite" cluster of filled pores, 2 is the pore from which the liquid is 
extruded, and 3 is the formed cluster of filled pores extruded from which is impossible because 
of the break of connection with the "infinite" cluster of filled pores. 

We will analyze successively all reasons for entrapment. To qualitatively analyze the energy 
reason for entrapment, we calculate a change in the energy of the system at the disappearance of 
a liquid nanocluster in one of the pores in the completely filled medium at Q\ = 1, which consists 
of pores of the same radius, when the overpressure vanishes (Figs. [2^, |2b). Following [6j, we 
accept that the chemical potential of the liquid in pores with a size of 13 nm does not change. 
In this case, a change in the energy at the extrusion of the liquid from a pore is determined by 
the difference between the surface energies after (E Q ) and before (Ei) the extrusion of the liquid. 
Energy Ei is the energy of the liquid-solid interface in the pore with throats the number of 
which is equal to the number z of the nearest neighbors: E^ = <Ji s (S — zS\), where 5* is the area 
of the pore, Si is the area of the throats, and a\ s is the solid-liquid surface tension coefficient. 
Energy E is the sum of the energy of the solid-gas interface a sg (S — zS\) (a sg is the solid-gas 
surface tension coefficient) and the energy of menisci formed in the throats of the pore <7\ g zS\ 
(crig is the liquid-gas surface tension coefficient, Figs. [2^, |2b). Thus, a change in the energy of 
the system at the extrusion of the liquid from the pore in the completely filled porous body 
AE(6i = 1) is given by the expression 

AE(6 1 = 1) = OigzSt + (a sg - a ls )(S - zSt). (1) 
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Figure 2: Schematic representation of a change in the state of the porous medium and the 
number of menisci in the pore before and after its extrusion of the liquid from it at the filling 
factors of the porous medium (a, b) 9\ = 1 and (c, d) 9\ < 1 and (d, e) schematic representation 
of the cluster of filled pores in the absence of paths for the extrusion of the liquid 

For the liquid that does not wet this porous body, we have a sg < cr/ s ; for this reason, the 
extrusion of the liquid from the pore is energetically profitable if AE(9 = 1) = a\ g zS\ + (a sg — 
&is)(S — zSi) < 0. Correspondingly, if a\ g zS\ > (ai s — <J sg )(S — zSi), the liquid is not extruded 
from the pore in the completely filled porous body. Estimates of AE show that entrapment is 
possible because the energy spent on the creation of menisci in the throats of the neighboring 
pores is higher than the energy released at the disappearance of the porous body-nonwetting 
liquid interface. 
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In a partially filled porous medium, 9\ < 1 and a filled pore in the initial state is surrounded 
by both the filled and empty pores (Fig. [2b) with menisci in their throats. When the liquid is 
extruded from such a pore, these menisci disappear, but new menisci are formed (Fig. [2U). If 
the pore in the initial state has k menisci, the number of menisci in the final state is z — k. As 
a result, a change in the energy of the system at the extrusion of the liquid from the pore in 
the partially filled porous medium is 



Therefore, if <J\ g {z — 2k) Si > (cr^ — <J sg )(S — zSi), the liquid nanocluster is stable and the 
liquid is not extruded from the pore in the partially filled porous body. 

According to Eq. (refeq2), the condition of entrapment depends on the number of empty 
pores k around the filled pore. A randomly chosen pore in the partially filled porous medium 
is in contact with k empty pores and, therefore, is surrounded by z — k filled pores. We assume 
that « 1 on average; in this case, k « (1 — Q\)z. The condition of entrapment of the liquid 
in the pore becomes 



where 77 is the connectivity factor, which is the ratio of the area of menisci average in the 
environment to the pore surface area S, and S m is the area of all menisci of the pore. In the 
general case, S m is the area of the menisci per pore in all throats of surrounding pores. 

Under the accepted assumptions, inequality ([3]) is not satisfied at 9\ < 1/2. Consequently, 
entrapment is possible only when the filling factor of the porous medium is Q\ > 6 c q ~ 1/2. At 
8\ < $c(b the number of menisci in the initial state is larger than that in the final state, and 
the process in which these menisci disappear and the number of menisci newly appearing in 
the neighboring pores at the extrusion of the liquid is smaller than the number of disappearing 
menisci is energetically profitable. Thus, the state of the (nonwetting liquid-nanoporous body) 
system is metastable at 9\ < 8 c o; i.e., after removal of the overpressure, the extrusion of the 
liquid from the porous medium is energetically profitable. At 6\ > 8 c q, the stable state of the 
liquid is a state in which stable liquid nanoclusters are present in pores. In the general case, 
8 c o is determined by the structure of the porous medium, in particular, depends on the pore 
size distribution. 

A change in the energy of the system at the extrusion of the liquid from a randomly chosen 
filled pore (the disappearance of the liquid nanocluster in the porous medium) after removal 
of the overpressure can be calculated using analytical percolation theory methods [23] ■ To this 
end, it is necessary to calculate a change in the energy of menisci average over all possible 
configurations (mutual arrangements) of pores of various sizes. We perform this calculation in 
the mean field approximation for the geometrical state of the medium with the infinite cluster of 
filled pores. A filled pore with a liquid nanocluster has neighbors, and contacts (interactions) 
between nanoclusters occur in the throats of pores, where a meniscus appears. Energy AE 
that should be spent on the extrusion of the liquid from the pore with radius R in the medium 
with filled pores of radii {Ri} at pressure p depends on the filling factor of the porous medium 
and radii R and {Ri}. This energy is the sum of the energy of the porous medium-liquid 
interface Sei(R, {Ri}) and the energy necessary for the formation of menisci in the throats of 
the neighboring pores 5ei nt (R, {Ri}, #1), which is due to the interaction of the liquid nanocluster 
in the pore with liquid clusters in surrounding pores: 




(2) 




(3) 
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AE(R, {R,}, 9,) = 5e x {R, {R t }) + 6e int (R, {R t }, 9 l ) 



5e 1 (R,{R*})= P V-5a(l-r ] (R,{R t })S, r, = M^lIM, (4) 

5s int (R, {R^, 9,) = a5S m (R, {Ri}, Oi) 

Here, V is the volume of the pore, a is the surface energy of the liquid, 5a = (<Ji s — a sg ) is 
the difference between the surface energies of the solid-liquid and solid-gas interfaces, 5S m is 
a change in the surface of menisci at the intrusion of the liquid into a pore with radius R, and 
i] is the connectivity factor. Below, we will use the mean field approximation; in this case, 

Se.iR) = P v- Mi - v(R, v(R, *i) = igrs^gi^iW)) , 

S (5) 

5e int (R,8i) = °(5S m (R, {Ri}, Oi)) 

Here, ( ) means averaging over all pores of the medium except for the given pore with 
radius R. In this approximation, the area of the menisci S m and its change SS m are determined 
by the nearest neighbors of the given pore. Therefore, 



1 f°° 

v(R, 0i) = J o z(R, Ru Oi)s m (R, RjfiRJdR! 
{6S m (R,{R i },6 1 )) = J dR 1 f(R 1 )z(R,R 1 )s m (R,R 1 )W 1 (z(R,R 1 ),9). 



(6) 



Here, z(R,Ri), s m (R,Ri) are the number of the nearest neighbors and the area of the 
meniscus of the given pore with radius R, respectively, and Wi(z(R, R\), 9\) takes into account 
the interaction between liquid nanoclusters in neighboring pores and is the difference between 
the average numbers of menisci before and after the extrusion of the liquid from the pore 
per nearest neighbor. The product of Wi(z(R, Ri), 9\) and the surface energy of the liquid 
in menisci determines a possible change in the energy of the interaction of liquid nanoclusters 
with the environment at the transition to the unstable state. 

The extrusion of the liquid from the filled medium first occurs at a decrease in the pressure 
through the formation of individual empty pores and clusters of empty pores near the percola- 
tion threshold on the shell of the infinite cluster of filled pores. At 9\ > 9 C , Wi(9\ > 9 C ) should 
be determined as the difference between the numbers of menisci after and before the extrusion 
of the liquid from a pore on the shell of the system of empty pores: 

W 1 (R,R 1 ,9 1 > 9 C ) = ^(1 - 9 l ) n {P{9 l )y- n Z ~ 2n - 

^ z ri\(z-n)\ (7) 

z = z{R,R 1 ,9 1 ). 

Here,P(#i) is the probability that the filled pore belongs to the infinite cluster of filled 
pores. The first factor is the probability that an empty pore is located near the infinite cluster 
of filled pores. The second factor is the probability that this pore is surrounded by z — n empty 
pores and, hence, has z — n throats. The third factor is the difference between the relative 
numbers of menisci after {z — n) and before (n) the intrusion of the liquid into the pore. The 
fourth, combinatory, factor takes into account the variants of the distribution of n menisci in 
the nearest neighbors of this pore. Summation takes into account all possible variants of the 
mutual arrangement of empty and filled pores. 
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It follows from Eq. (J7|) that Wi(#i) at Q\ — Y 1 tends to 1, which corresponds to a change 
in the number of menisci at the extrusion of the liquid from one pore in the completely filled 
porous medium. A decrease in the filling factor is accompanied by a decrease in W\{Q\) and 
a change in sign at Q\ = 0o(z), which corresponds, according to Eqs. and (jH]), to a change 
in the sign of the energy of the interaction of liquid nanoclusters with the environment at the 
transition to an unstable state. In this case, liquid nanoclusters in pores can become unstable 
and be extruded. At 9\ < 6 C , Wi(9i) is constant, which corresponds to the destruction of 
individual liquid nanoclusters. 

The connectivity factor rj(R, 9{) and quantity (5S m ) can be calculated using a particular 
model of the porous medium. In this work, they are calculated using the model of randomly 
located spheres j20] with pores of various radii taking into account correlations in the spatial 
arrangement of pores in the medium. In this model, 

I r\R+Ri\ 

z(R,R 1 ) = —— g 1 (R,R 1 ,r)-d 3 r, (8) 

where V pore is the volume of one pore and 

gi{R,R u r) = <p H i 

(9) 



R 2 -(r- R 1 ) 2 n n 

x— , y = R + R\—x — r 

2r 

is the pair correlation function of the system of randomly located spherical pores with radii 
R and and 

\2\ 



1 r R+R i nR (R\ - (r - Rf) 

VJ\R-K\ r (1Q) 

V = ^-((R + R 1 ) 3 -\R-R 1 \ 3 ). 

According to calculations by Eqs. ()6]) and (18|)-( IT0|) . 77 ~ (Ro/R)~ a ,a ~ 0.3 in this model, 
where Rq is determined by the pore size distribution function and is about the minimum radius 
of pores. In this case, since the pore size distribution is narrow, 

(8S m (R,{Ri},9)) = 

= J dR 1 f(R 1 )z(R,R 1 )s m (R,Ri)W 1 (z(R,R 1 ),6)^ V (R)W 1 (z,e), 

where z is the average number of the nearest neighbors in the model of randomly located 
spheres, which depends only on porosity: z = — 81n(l — (p) |20] . 

It follows from Eqs. (EJ) and (jS]) that a change in the energy of the pore AE(9, R) = 
5ei(R,6i) + 56i nt (R,6i) after the extrusion of the liquid from it in the partially filled porous 
body depends on the radius of the pore R and the filling factor Q\. 

The possibility of the existence of liquid nanoclusters in pores in the stable state after 
removal of the overpressure is determined by the fluctuation probability w of entrapment of 
the liquid in the pore [16J. Since a "quantum" of quasistatic filling is filling of one pore and, 
therefore, a pore can be in one of two possible states, filled or empty, we have 

AE(9, Tj(R), R) ! 
w = (l + exp( )) . (11) 
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If AE (8,rj(R), R) > 0, the probability is w ~ 1 and the liquid nanocluster in the pore is 
stable, whereas if AE (9, rj(R),R) < 0, the probability is w = and the liquid is extruded from 
the pore. At the filling factor 9 > 9 c o, entrapment of the liquid is possible in pores with radii 
smaller than the critical value R*(8i,T) = Rq(1 + j^W^z, 9i)) «, which is determined from 
the equation AE = and depends on the filling factor 9\ and temperature T. Temperature 
dependence appears owing to the temperature dependences of the surface tension coefficients 
(Tig, ais, <j sg . The volume AV(9i) of the liquid confined in the porous medium filled to the filling 
factor 9\ is determined by the integral of the distribution function f(R): 

pR*{B x> T) 

AV{9 X ) = / R 3 f(R)dR. (12) 

Jo 

Further, we calculate the volume of the remaining liquid determined by the pore space 
geometry. At the filling factor 9 < 1, the probability W v (9i) that the pore belongs to a 
geometrical cluster of filled pores that is not connected with the infinite cluster of such pores is 
determined by the product of the probability 9i — P{9i) that the pore does not belong to the 
infinite cluster of filled pores and the probability 1 — 9% that the surrounding pores are empty. 
This condition guarantees the absence of paths for the extrusion of the liquid. If the number 
of the nearest neighbors in the porous medium is z, this probability is given by the following 
expression, where all possible geometrical configurations of clusters of n filled pores are taken 
into account: 

= $> - pwni - e 1 y- zl = (i-fli)'-a-mr (13) 

' zn\\z — n ! z 

n=l v ' 

Here, all possible configurations of the arrangement of filled and empty neighboring pores 
are taken into account. The total volume of the liquid remaining in these clusters is determined 
by the fraction of the pores ANi(9\) located in the clusters in which the liquid remains because 
possible paths for extrusion are absent: AiVi = Jg 1 d9W v (9). According to Eq. (TT3|) . ANi(9i) 
is independent of the temperature and, above the percolation threshold (when 9\ > 9 C and 
P{9\) > 0), tends to the constant AiVoo = L d9W v (9), which is determined by the dependence 
P(9i) and the number of the nearest neighbors z. 

Finally, the total volume of the liquid remaining in the porous body filled to the filling 
factor 9\ is given by the expression 

rR*{0i,T) roc 

AV(9 1 ,T)= / R 3 f(R)dR + / AN 1 (9 1 )f(R)R 3 dR. (14) 

Jo Jo 

Using Eq. (|T4|) . we performed calculations for the description of the experimental data 
obtained. The lines in Fig. [1] are the dependences of the fraction of the nonwetting liquid 
confined in the porous body on the filling factor at various temperatures as calculated by 
Eq. (Il4p and the temperature dependence of the fraction of the nonwetting liquid remaining 
in the porous body at the filling factor 9\ = 1. The calculations were performed with the 
parameters of the porous medium and liquid (<p, R,5R/R) corresponding to the experimental 
values from [20J . The surface tension coefficient of water and its temperature dependence were 
taken from [24J. The surface tension coefficient of water at T = 293 K is 72 mJ/m 2 [25]. The 
quantity 5a and its temperature dependence were determined from the temperature dependence 
of the extrusion pressure using the method described in At T = 293 K, 8a = 10 mJ/m 2 . 

As can be seen in the figures, the resulting theoretical dependences well describe the entire 
set of the experimental data, which can indicate that the proposed model is applicable not 
only to the description of entrapment in the system studied in [19], but also to the description 
of the intrusion-extrusion process in other systems for revealing the structural features of a 
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porous medium. The applicability of this model to systems of various natures and for various 
experimental conditions will soon be analyzed. 

Further, we discuss the stability of the state of the nonwetting liquid in the porous body. To 
this end, we calculate the addition to the thermodynamic potential of the (nonwetting liquid- 
nanoporous medium) system owing to the dispersion of the liquid. This addition is obtained 
by averaging the energy AE(8, R) = Sei(R, 9\) + 5ei n t(R, 61) over the sizes of pores taking into 
account the probability of entrapment of the liquid in a pore with radius R given by Eq. ([8]) at 
the filling factor 81 and by multiplying the resulting expression by the fraction of pores in the 
porous medium: 



Estimates of AF(8i,T) at temperatures T = 279, 286, and 293 K show that the induced 
addition to the free surface energy of the system of liquid nanoclusters owing to the dispersion of 
the nonwetting liquid is small and, for the L23-water system, is AF(6i = 1, T = 293) ~ 3 • 10 -2 
eV per pore at temperatures above T c f» 287 K, which corresponds to the total extrusion of 
the liquid from the porous body for the L23C8-water system under consideration. At lower 
temperatures, AF{9 1 = 1,T = 286) « 0.3 eV and AF(0i = 1,T = 279) « 1.6 eV. Thus, the 
nonwetting liquid in the nanoporous body at T < T c can pass to a stable dispersed state; i.e., it 
efficiently becomes "wetting". Extrusion is energetically unprofitable for those clusters of filled 
pores for which the number of menisci formed in the final state after the extrusion of the liquid 
is larger than the initial number of menisci. In this case, the extrusion of the liquid should be 
accompanied by an increase in the surface energy of the liquid. For this reason, entrapment 
can be treated as the attraction between liquid nanoclusters in the neighboring pores. A 
detailed analysis and comparison of the proposed model with all reported experimental data 
on entrapment of the liquid in porous media will be given elsewhere. 
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